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1. INTRODUCTION 


The basic reasoning behind the epidemiological 
studies of any disease is to give guidelines to plan- 
ners about the importance of the disease in an area 
and to monitor and evaluate changes in the en- 
demicity due to control programmes. Since most 
of the disease control operations are long term 
programme and implemented in a phased man- 
ner, information in every phase of the programme 
(whether it is a planning, preparatory phase, 
evaluation or implementation phase) need to be 
processed, quantified and rational feedbacks 
need to be given for the improvement if any re- 
quired. The disease could be quantified by es- 
timating different parameters. These parameters 
can be estimated either by direct crude methods 
or by use of sophisticated mathematical formula- 
tions. Estimates obtained by either by crude 
method or by mathematical model can be used to 
evaluate the effectiveness of control programmes 
during different phases. Therefore, adequate 
knowledge on statistical/mathematical principles 
of designing studies, collection of correct and ade- 
quate sample and a rational analysis, etc. are high- 
ly essential not only for carrying out the study but 
also for drawing meaningful conclusions. 


Though the theory of mathematics developed for 
biological applications are extensive and often so 
abstruse, its implications are no longer accessible 
to many biologists. Therefore, it is an essential 
task for applied mathematicians to maintain a 
communication link between the mathematical 
theory and biological reality. Such a communica- 
tion link, if rightly understood, is expected to 
result in far reaching practical consequences. In 
this training module the authors tried to develop 
a link between mathematical theory and biologi- 
cal reality in relation to an important vector borne 
disease filariasis. 


14. FILARIASIS (EXTENT OF THE 
PROBLEM): 


Filariasis is a major health problem in tropical 
countries of the world. Though mortality due to 
filariasis is negligible, there is a high degree of 


morbidity due to its acute and chronic manifesta- 
tions. According to latest WHO estimates 2677 
million people are living in countries where the 
disease is endemic, 905 million people are at the 
risk of infection (living in areas where transmis- 
sion is active) and 90.2 million people are in- 
fected'*, Over two thirds of lymphatic filariasis 
cases of the world are contributed by India, China 
and Indonesia. Lymphatic filariasis is an age old 
disease recurded in India as early as 6h century 
B.C. Though both bancroftian and brugian forms 
of lymphatic filariasis are prevalent in India, 
bancroftian is found to be the commonest, ac- 
counting for 98% of the filarial infection. A con- 
tinuously increasing trend was observed from 
1962 (5.3 million mf carriers) to 1989 (27 million 
mf carriers). Published results have shown that 
377 million people in India are exposed to the risk 
of W.bancrofti and B.malayi infection with an es- 
timated 27 million mf carriers and 19 million 
chronic filariasis cases. 


The study of transmission dynamics of filariasis is 
complex involving man, vector and parasite and 
several attempts have been made to study the 
relationship between host, vector and parasite. 
However due to multiplicity of variables and 
processes which are involved makes it difficult to 
arrive at a comprehensive quantitative under- 
standing of the epidemiology and control of 
filariasis without the aid of mathematical/ statisti- 
cal tools. 


In the present communication we describe the 
statistical methods useful for collecting, analyzing 
and interpreting data for understanding 
epidemiology of filariasis and its control. The 
methods that are currently used have been 
published in the form of manuals and books , 
In this manual we confine ourselves to the new 
mathematical/statistical methods that are 
developed for analyzing the large data base at the 
Vector Control Research Centre, Pondicherry. 
This is expected to meet the needs of not only 
planners and implementor of national filariasis 
control programmes but also researchers engaged 
in filariasis studies. It is hoped to stimulate both 


biologists and applied mathematicians/ statis- 
ticians to improve our present level of under- 


standing of the dynamics of lymphatic filariasis. 


2. COLLECTION OF PARASITOLOGICAL DATA 


21. THE NEED FOR A SAMPLING 
DESIGN: 


In any epidemiological investigation it is necessary 
to measure the prevalence or incidence of the dis- 
_ ease. Theoretically this could be done by examin- 
ing all the individuals in a population. Such 
measurements are generally called ‘census’. A 
complete survey of all the individuals in a com- 
munity will be useful not only to detect and treat 
the carriers but also to have an idea of the extent 
of the problem under study. As for as the first ob- 
jective is concerned, it will be fulfilled, provided 
at least a majority of the population is sampled, if 
not a 100 percent survey. However 100 percent 
coverage is impracticable when the number of in- 
dividuals in the population is very large. The com- 
- plete survey of individuals for meeting the second 
objective is ‘of limited value in terms of not only 
both cost and time but of less precision of the es- 
timates. In such cases a small collection from the 
larger population will be studied and inference 
can be drawn about the population. This process 
is usually called ‘sampling’. A survey following a 
properly designed probability sampling can 
provide reliable estimates of disease prevalence 
or incidence with a greater precision within the 
available resources. : 


The conventional blood surveys, in filariasis con- 
trol programmes, are neither census nor sample 
survey. These surveys are aimed at collecting as 
many slides as possible, without resorting to any 
specific method of data collection. The main ob- 
jective of mass blood survey is to identify all the 
parasite carriers in a community and to eliminate 
the foci of infection either by drug therapy or by 
isolating the carriers. Using such data for measur- 
ing prevalence or incidence of microfilaria in the 
population may give erroneous picture due to its 
non-random sampling procedures. These surveys 
are subject to considerable amount of errors 


which could be avoided by adopting the scientifi- 
cally more valid sampling techniques. Further, it 
does not provide any objective basis for measur- 
ing the precision of the sample results. The 
precision is determined and measured in terms of 
the standard deviation of the sample estimate. 
The commonest error is the indiscriminate use of 
‘nqin’ (where ‘p’ is the proportion of positives, 
q =1-p and ‘n’ is the sample size) as the variance 
of a proportion as a measure of precision, without 
regard to the manner in which the sample was 
selected. If the distribution of disease is very 
sporadic or concentrated in a small part of the 
region, a sampling design based on the prob- 
ability distribution of the disease would be more 
suitable than a mass blood survey based on 
stratified random sampling of individuals. 


Census or total survey is neither feasible nor 
necessary for estimating disease parameters, 
mainly because of cost and time. The more reli- 
able and scientific estimates of the parameters 
based on sample survey can considerably reduce 
the cost and time. However, the adequacy of the 


- ‘sample and reliability of estimates derived from it 


are guaranteed only if carefully selected prob- 


ability sampling designs are used. Although many 


sample designs are available for a particular 
problem, it is necessary to make a rational choice 


among them. A rational choice can be made only 


if comparative efficiency of various sampling 
designs are known. Eventhough the design of 
sampling may vary from place to place depending 


- onthe objective of the study, as ageneral guideline 


the following steps must be followed to get reliable 
BIEps. ; 


2.1.1. Principal steps involved in a sample survey: 
a) Define the objective of the survey: 


Since the sampling design and the entire 


procedure rests on the objective of the survey it is 
important to define the purpose of the survey (e.g. 
Estimate the prevalence/ incidence or detect all 
parasite carriers etc.). 


b) Define the population to be sampled: 


The population to be sampled should be clearly 
defined and it depends on the objectives of the 
study. For example if the objective is to evaluate 
the effect of vector control on filariasis \ransmis- 
sion, then a longitudinal follow up of the popula- 
tion of children between 0-8 years of age are the 
population to be sampled. However, if the objec- 
tive is to identify all the parasite carriers for 
chemotherapeutic control, then a complete survey 
of the population is necessary. If the study is of 
analytical in nature (eg. studying the prevalence, 
to monitor and evaluate impact of control 
programmes in the population) then a cross sec- 
tion of the entire population can be sampled. 


c) Data to be collected: 


Proper care should be taken to include all the 
relevant information for the study. Relevant in- 
formations are those that may have impact on the 
study variables either directly or indirectly. 


_ For example: ‘An epidemiological evaluation of 


the effect of vector control programme for 
filariasis can be measured by assessing the infec- 
tion prevalence and intensity of microfilaria 
before and after control. Since prevalence of in- 
fection in a locality depends on many factors such 
as age and sex distribution of the population, 
mosquitogenic conditions in the locality, 
socioeconomic conditions, family size and method 
of control etc., it is essential to collect information 
on the above factors. Collecting unwanted and ir- 
relevant information may lead to unnecessary 
confusion and thus affect the quality of the study. 


d). Degree of precision desired: 


Since only a part of the population is surveyed in 
sample study, the results are subjected to some de- 
gree of uncertainty. Though it can be reduced con- 
siderably by taking large samples and by using 
more sensitive methods of measurement (eg. 


membrane filter technique, ELISA techniques for 
parasitological surveys in filariasis studies), it 
needs large amount of money and takes more 
time. Consequently, based on the degree of 
precision required in the results, sample size will 
vary. This can be fixed by the researcher in con- 
sultation with a statistician. 


¢) Methods of measurement: 


Several methods may be available for measuring 
the variables under study. The researcher should 
be well aware of the limitation of the different 
techniques available for measurement. 


For example, a filariasis parasitological survey can 
be carried out by finger pricking (either 20 cmm 
or 60 cmm), membrane filter technique and 
ELISA. Membrane filter technique and ELISA 
may be more sensitive techniques but in large 
scale surveys they are less cost effective and time 
consuming. On the other hand finger pricking is 
both cost effective and less time consuming, it is 
not a sensitive method to detect very low 
microfilaria carriers. Hence a standard method 
of collecting information is always advisable for 
any study. A major part of the preliminary work is 
the construction of record forms on which the 
questions and answers are to be entered. 


f) The frame: 


Defining the sampling unii is important before 
selecting the sample. These are the parts of the 
population which should not overlap with each 
other. These may be an individual or household 
etc. The construction of such list is called sampling 
frame, is one of the major practical problems in 
sampling. 


Example:In parasitological surveys each in- 
dividual is a sampling unit. Since construction of 
a sampling frame in which information about each 
individual is very difficult to obtain before con- 
ducting a survey, each family may be considered 
as a sampling unit. Hence it is sufficient to have a 
sampling frame consisting of information about 
the head of each family. 


g) The pretest: 


It is always advisable to test the questionnaire in 
the field on a small scale. This will give an idea 
about the cost and time required to complete the 
study. The problems encountered in this can be 
rectified to a greater extent in the main study and 
thus reducing the measurement bias. 


h) Organization of the field work: 


Before starting the survey the field staff should be 
trained with the objective and with the methods of 
measurement. Proper supervision of the field staff 
and their work will guarantee the collected infor- 
mation. Steps should be taken in advance to min- 
imize the non-response bias and to tackle such 
situations. 


i). Summary and analysis of data: 


The significant part of any study is the analysis of 
data. Proper application of statistics for obtaining 
results is important. Improper application of 
statistics may lead to wrong conclusion about the 
population. The following are the sequence of 
statistical operations to be carried out on any field 
collected raw data’: 


(a) Investigation in to validity of assumptions, 
biases, etc. must be carried out. 


(b) Descriptive statistics: The descriptive statis- 
tics such as mean, standard deviation, can be 
presented in tabular form. Graphical repre- 
sentations can be included wherever neces- 
sary. 


(c) Standard indices: Indices such as infection 
prevalence, intensity, resting mosquito den- 
sity, etc. can be computed. 


(d) Parameter estimation: Some sensitive 
parameters can be estimated more accurate- 
ly with help of mathematical models. 


(e) Trends: Trend analysis may be carried out to 
see the changes in the parameters over a 


period of time (e.g., change in prevalence of 
microfilaraemia). 


(f) Correlation and Regression: Correlation 
analysis can be carried out to see the degree 
of relationship between dependent and inde- 
pendent variables and the regression analysis 
will be useful to identify the predictor vari- 
ables and to establish relationship, if any, be- 
tween the study and predictor variables. 


(g) Contigency: Contigency tables can be 
prepared to analyse the qualitative data using 
chisquare tests, etc. 


j) Information gained for future surveys: 


Any completed sample survey is potentially a 
guide to improve further sampling. The difficul- 
ties encountered, the level of error and other in- 
formation about the population will improve the 
precision of estimates in the future sampling. 


2.1.2. Use of supplementary information: 


Supplementary information may be used to 
provide greater precision in our sample estimates. 
Stratification, a method for utilizing supplemental 
information, based on past data may be used to 
divide the population into groups such that ele- 
ments in each group are more alike than are the 
elements in the population as a whole. In case of 
filariasis control programmes there exist a greater 
heterogeneity in the transmission potential of dif- 
ferent localities in an urban situation. This 
heterogeneity is the result of unplanned urbaniza- 
tion and hence some areas are more prone to the 
risk of transmission’, Any sampling procedure 
for filariasis without giving serious thoughts to the 
following factors would lead to over/under estima- 
tion of the parameters under study: 


a) The distribution of microfilaria is known to be 
overdispersed in all age groups of the human 
population. 


b) The microfilaria carriers are over crowded 
due to the focal transmission of the parasite. 


¢) The conventional blood sampling technique 
known to provide bias in the estimation of 
parameters of filariasis control programmes. 


2.1.3. Probability sampling methods: 


The conventional probability sampling methods 
are available in many text books®”'!"3. Here we 
describe those methods with their advantages and 
disadvantages while conducting surveys. A 
population can be sampled by using any one or 
combination of the following probability sampling 
procedures:. 


a) Simple Random Sampling (SRS) 
b) Systematic sampling 

c) Stratified random sampling 

d) Cluster sampling 


Simple Random Sampling: This is a method in 
which each individual in the population has an 
equal chance of being included in the sample. It 
needs a list of serially numbered identification of 
individuals in the population. Individuals are in- 
cluded in the sample by selecting a sample of num- 
bers from random number tables. SRS is simple 
in theory but is of less practicable when large scale 
surveys are undertaken. It needs a complete sam- 
pling frame of all the individuals in a population. 


Disadvantages of Simple Random Sampling 
(SRS) of individuals: 


1. Since in a SRS procedure, individuals are con- 
sidered as sampling units, the preparation of 
sampling frame which needs the information 
about each sampling unit is difficult when large 
scale surveys are undertaken. 


2. In SRS only selected individuals in a house or 
locality are surveyed. This method, which is 
restricted only to specific individuals in a 
house, creates suspicion among the individuals 
in a family which might affect the precision of 
the estimate due to the non response bias. 


3. Since filarial parasite carriers are known to be 
aggregated, the SRS of individuals as sampling 


units fails to give an estimate of aggregation of 
parasite carriers in the community. 


Systematic sampling: In this procedure the 
selected sampling units are spaced regularly 
throughout the sampling frame. A starting num- 
ber ’r’ is chosen randomly from a sampling frame 
and thereafter every r'” individual is selected. For 
example to get a sample of 500 persons from a 
population of 10000 we proceed as follows: 


(a) calculate the random interval, 
k = N/n = 10000/500 = 20 

(b) Select a number ‘r’ from 1 to k using random 
number tables. 

(c) Select every r'" unit from the sampling frame. 

(d) Continue step (c) until the required sample 
size 1s reached. 


Systematic sampling has three major advantages: 


i) The investigator does not need to know the 
sampling frame in advance; the frame can be 
constructed as the study progresses. 


li) It is simple to implement under field condi- 
tions. For example it iseasierto take blood 
smears from every fifth house in a block than 
to determine which houses are to be visited by 
means of random number tables. 


iii) If a trend is present in the sampling frame, 
such as from small units to large units, a sys- 
tematic sample will ensure coverage of the 
spectrum of units. 


For example: A systematic sample of hospital ad- 
missions selected over a period of time may 
provide a better estimate of prevalence of a com- 
municable disease than an estimate based on SRS. 


A potential disadvantage of using systematic 
sample is that it may provide a biased sample if the 
population contains a periodic type of variation, 
and if the interval between successive units in the 
systematic sample happens to coincide with the 
wave length. 


Stratified random sampling: In this procedure 
the population is divided into strata, or groups of 
units having certain characteristics in common, 
and a sample of units is drawn from each stratum. 
The samples from each stratum can be selected 
using any one of the probability sampling 
methods. The method of selection of samples 
need not be identical in all the strata. For example 
the investigator can apply the SRS technique in 
one stratum, cluster sampling in the next stratum 
and systematic sampling in another stratum etc. 
Accordingly the sampling procedure is called as 
Stratified SRS, Stratified cluster and Stratified 
systematic sampling procedure etc. Stratified 
sampling is widely used and has several ad- 
vantages over SRS. 


i) The investigator can ensure that each sub- 
group of the population is represented. 


li) Since the population is divided into 
homogeneous subgroups, more precise es- 
timates of the population parameters can be 
obtained. 


The main potential disadvantage of stratified sam- 
pling is that loss of precision can occur if very small 
numbers of units are sampled within individual 
strata. 


Cluster sampling: The individuals in population 
are often grouped into clusters (e.g.,families, vil- 
lages, hospital wards etc.). In such situations the 
cluster sampling is more suitable than any other 
sampling procedures. In this procedure a sample 
of clusters (e.g.,city blocks) are randomly selected 
from the population and then observations are 
made on all individual units (households) in each 
chosen cluster. Though sampling variance is 
greater for cluster sampling, it could be reduced 
by increasing the number of clusters to be 
sampled. This procedure has two major ad- 
vantages: 


i) Sampling frame is not required for the entire 
population. 


ii) Economic and less time consuming. 


Example: Suppose the objective is to study the 


prevalence of malaria or filariasis in a district. As- 
sume that villages are grouped to form clusters 
which are identical in all ecological parameters. 
Then a sample of clusters may be drawn and all 
the individuals in the chosen clusters are ex- 
amined. 


2.1.4. The sample design for filariasis survey: 


Considering the above factors it would be ideal to 
stratify the study area according to the 
tmicrofilaria prevalence or entomological observa- 
tions (zones). In the absence of parasitological 
and entomological information about the study 
area, stratification based on sanitary and living 
condition of the people may be considered. Asa 
second stage, each stratum may be further clas- 
sified into substrata (eg: streets) consisting of 
homogeneous sampling units. The sample size for 
each substratum has to be determined either using 
prior estimates of the population characteristics if 
available or in proportion to the population of 
each substratum. The sampling unit may repre- 
sent the households which forms a cluster of in- 
dividuals within the family. Further sampling may 
be carried out by adopting any one of the follow- 
ing schemes: 


hen mupling th fin: 


Each stratum is considered as a cluster of 
homogeneous streets (substrata). A sample of 
clusters are selected from each stratum based on 
simple random sampling without replacement 
(SRSWOR) technique. Since in this procedure 
only the substrata are sampled and all the in- 
dividuals in the selected substratum are surveyed 
this is called a single stage cluster sampling within 
the stratum. : | 


heme 2.Sampling the | 


In this procedure the sampling is carried out at — 


each substratum level. Each zone (stratum) is fur- 
ther stratified into substrata, consisting of cluster 


of homogeneous households. Then sample of 


households for each substratum are selected 


either by SRSWOR or Systematic sampling — 


methods. The selected households in all substrata 


are surveyed and each household is considered as 
a cluster. 


2.1.5, The Estimation of Sample Size: 


The method of estimation of sample size varies ac- 
cording to the availability of past information 
about the population characteristics. As stated 
earlier if no prior estimates of the population 
characteristics are available, the sample size for 
each stratum has to be determined in proportion 
to the population of that stratum. 


[f scheme 1 is decided to adopt as the sampling 
protocol then a sample of substratum is selected. 
The number of clusters (streets) to be sampled in 
each stratum is based on the number of 
households to be surveyed in the zone and the 
number of houses in the cluster. This procedure is 
carried out as follows. Based on the average fami- 
ly size, the number of households to be surveyed 
for a zone is calculated from the zone sample size. 
Let’a’ be the number of households to be sampled 
in the zone. Select a street within the zone ran- 
domly and note the total number of houses in that 
street. Continue the sampling of streets until the 
number of households of the selected streets 
equals to ’a’. This gives the number of streets 
(clusters) to be sampled for the zone. Then survey 
all the households in the selected clusters. 


In scheme 2, the sample size of the zone is allotted 
to each substratum according to their population. 
The number of households are determined based 
on the average family size of the street. The units 
(households) are randomly selected using random 
number tables. Survey is carried out only for the 
selected houses in all streets. 


The principal steps involved in the choice of a 
sample size are well described in many textbooks” 

~ Initially, the sample size is estimated based 
on the total population of the study area. Then, 
this is distributed to different strata and substrata 
according to their population size. Since 
microfilaria count/carriers are known to be high- 
ly aggregated and this distribution is adequately 
described by the negative binomial probability 
law, the calculation of sample size must take into 
account of this factor. The sample size is deter- 


mined by the following factors: 


a) Mean number of mf count/carriers per per- 
son/house (m). 

b) The negative binomial dispersion parameter 
(k). 

c) Desired level of precision (r) as a percentage. 

d) Probability of achieving the desired level of 
precision (q). 


The sample size can be estimated by the following 
formula: 


n = (100tq)*(1/m + 1/k)/t? 0... (1) 
Where n = Sample size required. 


tq = Students’s t value for n-1 degrees of 
freedom at a % level of significance. 


m= Mean number of mf carriers/mf 
counts 
k = m7/(s2-m) veccssseeecccarsseeeeee (2) 
s?= variance of mf carriers/mf counts. 
r= desired level of error (%). 


For most practical purposes we can assume 
tq =2.0 (for 95% confidence limits). Then the 
formula (1) is simplified to, 


n = 40000 x (1/m+ 1/k)/t® ener (3) 


Example: The mean mf count and variance based 
on the 1981 filariasis survey in Pondicherry are 
given as 0.9125 and 31.9378 respectively. From 
this, the value of ’k’ was estimated by using equa- 
tion (2) as 0.0268. Based on these informations the 
sample size for a filariasis study in Pondicherry 
can be estimated using equation (3) as 15342. 


This sample size can be distributed to the entire 
study area which was divided into six zones 
(stratum) with homogeneous sanitary conditions, 
according to their population size. The number of 
households to be sampled is determined based on 
the average family size and the street population 
accordingly. The households to be surveyed are 
selected using random number tables and then all 
the members in the selected houses are surveyed. 


There are many situations in which the prior in- 
formation about the population characteristics 


may not be available. In such situations either a 
pilot study may be carried out to have an idea 
about the estimate of the parameters or a 
deliberate selection of about 10% of the popula- 
tion may be sampled. 


2.1.6. Parameter Estimation: 


Scheme 1: 


Let ’mijk? be the number of individuals in the k"” 
household of the es street (cluster) in i” zone. 


Then, 


a) the street wise mf prevalence (pij), which is the 
proportion of individuals positive for 
microfilaraemia in the j'" cluster unit of i‘ 
zone is given by 


Pij = Laijk/>Mijk X 100 


b) the Zone wise mf prevalence (pj), which is the 
ratio of number of mf carriers in the i zone 
to total examined in that zone is given by 


Pi = LLaijx/ZXmijx X 100 
jk jk 


c) Over all mf prevalence is given by 


Pp = LLLaijx / LL Xmijx X 100 
ijk ijk 


Total no.of +ves in all zones 
1.€. P = ---------------------------------------- X 100 


Total number surveyed in all zones 


Where, 
aijk 1s the number of mf carriers in the k'" 


household of j'” street in i” zone. 


Structurally all these are typical ratio estimates 
which are slightly biased, though the bias is seldom 
likely to be of practical importance. The estimated 
variance of the overall prevalence ’p’ is given by, 


(1-f){Zai7-2pLaimj + p* Y mj*} 


v(p) = 
nm‘ (n-1) 


where, 


m = >mj/n 

= Average number of individuals per cluster. 
a; = Number of mf carriers in the i'” zone. 
n = No.of clusters and f = n/N. 


Scheme 2: 


Let ’mijx’ be the number of individuals in the k"" 
household (cluster) of i street in i'" zone. Then, 


a) the street wise mf prevalence (pij), which is the 
proportion of individuals positive for 
microfilaraemia in the a street of i'" zone is 
given by 


Ppij = Laijk/2>Mijk X 100 


b) the Zone wise mf prevalence (pi), which is the 
ratio of number of mf carriers in the i zone 


to total examined in that zone is given by 


pi = LLaijx/Z2Xmijk X 100 
jk jk 


c) Over all mf prevalence is given by 


p = LLLaijk / LL Umijx X 100 
ijk ijk 


\ 


Total no.of +ves in all zones 
A ee X 100 


Total number surveyed in all zones 


The estimated variance of the overall prevalence 
’p’ is as described in scheme 1. 


2.2. SAMPLING DESIGN FOR CLINICAL 
SURVEY 


A clinical survey conducted along with blood sur- 
vey would provide valuable informations relating 
disease and microfilaraemia positivity. The main 
purpose of clinical survey is to study the 
prevalence of disease and to estimate the overall 
endemicity (includes mf + ves and disease + ves) 
which is an indication of the extent of the filarial 


problem in a community. Recent studies at 

VCRC have used the clinical data for estimating 
the true incidence by correcting the loss of infec- 
tion which includes individuals who are apparent- 
ly negative but clinically positive ( see section - 2.4 
). Further the progression of disease in relation 
to microfilaria status can be evaluated (see sec- 
tion - 2.5). The following clinical parameters are 
useful in the assessment of disease occurrence in 
a community’: 


(a) the incidence of acute manifestations such as 
adenolymphangitis, epididymo-orchitis, 
filarial fever, etc. and 


(b) the prevalence of above acute manifestations 
as well as chronic manifestations such as lym- 
phoedema (and its grades), hydrocele, 
chyluria, etc. 


The modified method of Kimura et q/'* canbe ap- 
plied for estimating the incidence of acute 
manifestations’. 


To generate data which are comparable, the clini- 
cal survey should take care of the following: 

(a) Defined criteria of disease?!©!”, All 
physicians should follow same criteria and 
should be trained specifically prior to actual 
survey, if necessary. 


(b) Should include both acute and chronic 
manifestations and actual physical examina- 
tion should be carried out. 


(c) Data of clinical history is useful but should be 
analysed separately. 


Though it would be ideal to examine clinically all 


persons covered during the blood surveys, it may 
not be practicable due to the following reasons: 


(a) Clinical examination of individuals is time 
consuming. 


(b) Availability of physicians for conducting the 
survey. 


Hence a subsample of blood smeared individuals 


may be examined clinically and whose history may 
be recorded. The actual sample size required can 
be estimated, if prior estimates are available or a 
pilot survey can be carried out to have an initial 
estimate of disease rate. If ’P’ is the disease 
prevalence and ’n’ the number samples to sur- 
veyed and ’r’ the percentage of error allowed in 
’P’. Then the sample size ’n’ can be estimated as 
follows: 


Here P, Q and r all are measured in percentages. 
The above formula assumes that the disease car- 
riers are binomially distributed in the population. 


However if the distribution is highly aggregated, 
then the required sample size may be underes- 
timated. In this situation the sample size is es- 
timated by assuming an aggregated probability 
distribution such as negative binomial, lognormal, 
Sichel etc. 


n = 40000 x (1/m + 1/k)/r’.......(3) 


Where, 
m = mean number of carriers per house. 
k = the negative binomial dispersion parameter. 
=m /(s’-m) 
s~” = sample variance. 


If the distribution is of random in nature then the 
sample size has to be estimated by using 
parameter of the Poisson distribution. Since the 
mean and variance are equal for the Poisson dis- 
tribution the formula for determining the sample 
size is very simple and is given as below: 


n = 40000 / mr’ 


where ‘m’ is the mean of the Poisson distribution. 


The frequency distribution of households accord- 
ing to number of diseased in each house, is given 
in the following table. Estimate the number of 
households to be sampled for estimating disease 
prevalence with a relative precision of 20%. 


Number of diseased Number of houses 
persons 


509 


Ww S 
— 
>_> © 


SOO AIAMNaPRWNKO 
meee DON Hh NN LY 


Total 661 


mean number of carriers per house (m) 


= 0.4191 
s*= 1.1862 
k = 0.2289 


The number of houses to be sampled (n) is ob- 
tained using equation (3) as 675.5. Hence a min- 
imum of 676 households can be examined 
clinically. A subsample of 676 households may be 
randomly selected (adopting any one of the sam- 
pling protocols discussed earlier under 
parasitological data collection) from the list of 
households who have parasitologically examined. 
This way a parallel parasitological and clinical in- 
formation can be obtained for the same 
household. 


The following indices will be useful: 


No. +ve for disease # 
(a) Disease Rate* = --------------------------- X 100. 


Number examined. 
# Confirmed by specific criteria. 


* For estimating the total disease rate both acute 
and chronic manifestations based on clinical find- 
ings (not based on history alone) should be in- 
cluded. 
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(b) Endemicity rate 


No. mf +ve (or) Disease +ve (or) 
both +ve 


ws sg Pa a tale ee X 100 


No. examined both for mf and disease 


(c) Incidence per 1000 population can also be es- 
timated by modifying the method given in 
section 2.4. for parasite incidence. It is ideal 
to consider age group 5-15 years for estimat- 
ing disease incidence. 


Further analyses on disease prevalence in 
relation to age, sex and microfilaria status can 
also be done, since age and gender crucially 
influence the disease status!®, Further this in- 
formation will enable to study the stability of 
the disease in a locality and also to observe 
changes on time scale ue to introduction of 
intervention measures ~. 


(d) 


2.3. FREQUENCY DISTRIBUTION OF 
MICROFILARIA COUNTS ~ 


2.3.1. Introduction: 


Certain degree of uncertainty is common in 
biological experiments partly due to the un- 
avoidable errors in the measurement processes 
and partly due to inherent variability in the 
biological processes”, The basic problem of 
Statistics in theory and practice is to estimate the 
population values based on the samples collected 
with minimum error. Though several estimates 
based on deterministic models are available, they 
fail to take in to account of the random variation 
in the data. The probability distributions are use- 
ful tools to study the inherent variation in the 
population process which cannot be directly 
measured. A trend in recent research on parasite 
ecology has been the recognition of the impor- 
tance of dispersion patterns to the population 
dynamics of host-parasite association”’, Es- 
timates based on the probability distributions are 
widely used in biology because these estimates are 
more realistic and gives better understanding of 
the population dynamics. Statistical models not 


only help to improve our understanding of 
population dynamics, but also provide a powerful 
tool for the assessment and evaluation of different 
approaches to control. If used sensibly, they can 
generate both qualitative and quantitative 
guidelines for the application purposes. In recent 
years fitting of suitable mathematical model which 
best fits the data finds application in most of the 
areas. The basic assumption in fitting of Statistical 
distributions is that the variable under study is ran- 
dom in nature with a definite probability of occur- 
rence. The extensive application of frequency 
distributions in biology is due to its accuracy and 
its reliability in prediction purposes. The quan- 
tification of sampling errors in parameter estima- 
tion which are common in biological phenomenon 
is the main feature of the probabilistic distribu- 
tions. Fitting a suitable distribution to the data 
does not convey any meaningful information un- 
less the rationale behind the probability distribu- 
tion is biologically justified. 


2.3.2. Frequency distribution: 


The term frequency is defined as the number of 
occurrence of an event in an experiment. This is a 
method of summarizing the field collected raw 
- data in the form of table. Such tables can be rep- 
resented by some known mathematical formulae 
which is unique for a particular biological or 
physical phenomena in nature. Usually the fre- 
quencies are expressed as the proportion of oc- 
currence to the total number of occurrences. This 
is also called relative frequency and total of all 
relative frequencies should be equal to one. In 
statistics the distribution of these relative frequen- 
cies are called the probability distributions. 


Distributions in nature are either continuous or 
discrete. Height, weight, wing length of an insect, 
instar duration etc are examples of continuous 
distributions. The normal, exponential, gamma, 
log-normal etc. are statistical distributions which 
are used to describe the continuous type of 
measurements. The discrete type of data such as 
number of eggs in an egg raft, microfilaria counts 
ina slide, number of infected mosquitoes, number 
of larvae in a dipper sample etc. can be described 
statistically by binomial, Poisson and negative bi- 
nomial which are discrete probability distributions. 
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Plants and animals in any geographical area are 
scattered about in one of three spatial patterns: 
Uniform, random or aggregated. Known mathe- 
matical frequency distributions canbe used as 
models for describing these spatial patterns””. If 
the frequency distribution of data from a sample 
fits one of these models, then: 


(a) the spatial dispersion of the population can 
be described in mathematical terms. 


(b) errors of population parameters can be es- 
timated. 


(c) temporal and spatial changes in density can 
be compared. 


(d) the effect of environmental factors can be as- 
sessed. 


The important distributions which are most com- 
monly used for describing discrete data are 
Binomial, Poisson and Negative binomial. An in- 
itial guess about the form of distribution is ob- 
tained by comparing the mean and variance (in 
other words the dispersion pattern of the data) of 
the observed frequencies. A simple measure of 
the degree of dispersion is the variance to mean 
ratio of parasite numbers per host”. The sample 
values are said to follow Binomial distribution if 
its mean value is significantly greater than the 
variance. If the mean and variance are ap- 
proximately equal then the data is supposed to fit 
Poisson distribution. There are certain situations 
where the variance is greater than the mean. In 
such situations the negative binomial distribution 
fits the data which is to represent the huge varia- 
tion in the datu. This is the most flexible of several 
possible distributions when the variance is sig- 
nificantly greater than the mean (o Se yw). 


2.3.3. Negative Binomial Distribution: 


Several mathematical distributions have been 
used as models for describing aggregated or 
clumped population and the negative binomial 
distribution is often a suitable model for inver- 
tebrate population. The negative binomial is the 
counterpart of the positive binomial and therefore 
the probability series of the negative binomial is 


given by the expansion of (q-p)* where p=,p/k 
and q=1+p. The parameters of the distribution 
are the arithmetic mean m, and the exponent k. 
The reciprocal of the exponent k (i.e 1/k) is a 
measure of the excess variance or clumping of the 
individuals in a population. As 1/k approaches 0 
and k approaches infinity, the distribution conver- 
ges to the Poisson series (o*—> ys ) Conversely, if 
clumping increases, 1/k approaches infinity 
(k— > 0) and the distribution converges to the 
logarithmic series. The general term of the nega- 
tive binomial distribution is , 


x+k-1 
P(x) = : (w/(w+k))* (k/(k+p))*.. (1) 
-1 


where, P(x) is the probability of observing x counts 
in a sampling unit. The k and mare the parameters 
to be estimated. Negative binomial can be applied 
to a wide variety of contagious distributions such 
as, 


a) True Contagion: The presence of one in- 
dividual or event increases the chance that 
another will occur in the same time. 


Constant birth-death-immigration rates: The 
growth of a population with constant rates of 
birth and death per individual and of immigra- 
tion per unit of time leads to a negative bino- 
mial for population size. 


b) 


c) mps: If clumps of in- 
dividuals are distributed at random and the 
number of individuals in the clumps are dis- 
tributed independently in a logarithmic dis- 


tribution, then a negative binomial will result. 


d) | mpling: If the 
population distribution is of continuous type 
and if the process of sampling is Poisson, then 
the resulting distribution is a compound Pois- 
son or negative binomial distribution. 


2.3.4. Method of Fitting Negative Binomial Dis- 
tribution: 


The fitting of negative binomial distribution invol- 
ves estimation of the parameters y and k where, 
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uy = the mean of the distribution and 
k = the dispersion parameter 

The parameter m is estimated by the mean’ m’ of 
the sampling distribution and an approximate 
value of ’k’ can be obtained by the following for- 
mula: 

k = m’/(s’-m).....(2) 
where m and s” are the mean and variance of the 
observed sampling distribution. An accurate es- 
timate of k can be obtained by substituting the 
value of ’k’ and’ m’ in the maximum likelihood 
equation. 


nin(1 + m/k) = ¥ (A@o/(k +x)) .... ..3) 
where A(x is the total number of counts exceed- 
ing x and n is the total number of sampling units. 
The probability values for various counts are cal- 
culated by substituting the values of m and k in 
equation (1). 


EXAMPLE: 
The following are the mf counts distributed in 


human population collected from an area in Pon- 
dicherry. 


0 300 
1 95 
2 70 
3 50 
4 15 
5 7 
6 - 
7 4 
8 3 
9 2 
10 1 
TOTAL 551 


In this example, 


1) since the distribution of mf within human is 
random, collecting blood smears from dif- 
ferent individuals form a mixture of Poisson 
process which leads to some aggregated dis- 
tribution. 


2) Mean is less than the Variance. 
mean = 1.0599 < variance = 2.5681 


Hence, the observed frequency distribution may 
be adequately described by the negative Binomial 
distribution which is the suitable distribution for 
representing the large amount of variation in data. 


The approximate value of k is obtained from equa- 
tion (2). 


ie k = 1.0599 *2/(2,5681-1.0599) 
= 0.7425. 


However, a real estimate of ’k’ can be obtained by 
trial & error method. i.e. by taking the values of k 
from 0.5 to 1.0 with increment 0.05. For each value 
of k the LHS and RHS of equation (3) can be es- 
timated. The difference between LHS and RHS 
can be obtained and tabulated as below. 


k m LHS RHS__ DIFFE- 
RENCE 
0.5 10599 62591 660.69  -33./9 
0.6 1.0599 560.69 568.73 -8.05 
0.6455 1.0599 535.31 535.77 -0.46 
0.649 1.0599 533.46 533.42 0.04 
0.8 1.0599 464.86 450.08 14.78 


LHS - Left Hand Side 
RHS - Right Hand Side 


Hence the value of k at which the sign of difference 
changes first time can be taken as the real value of 
i 

Then the estimated m and k values are substituted 
in the equation (1) to get the probability values for 
different counts. The expected frequencies are 
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calculated by multiplying the probability value 
with the total frequencies. A test of significance 
based on chi square distribution can be obtained 
by comparing the observed and expected frequen- 
cies. 


2.3.5. Frequency Distribution of Parasite in 
Human host: 


Helminth life cycles are often complex, involving 
more than one host species and many distinct 
Parasite development stages. The distribution of 
parasite numbers per host is important for under- 
standing the transmission dynamics of the parasite 
and also in the implementation of control 
policies, The frequency distribution of 
parasite numa per host have generated con- 
siderable interest in helminthology 23° The 
majority of distributions have been found to be ag- 
gregated or over dispersed and ecological theory 
strongly suggests that this serves to enhance the 
density dependant regulations of parasite abun- 
dance“! . Empirical evidences suggest that often 
the negative binomial probability distribution has 
proved to be a good peeenaie of observed pat- 
tern of parasite aggregation~ . Although, lognor- 
mal probability distribution can be used to model 
the observed pattern, it is less flexible to further 
asia treatment. Frequency distributions 
are primarily of value in the elucidation of qualita- 
tive as opposed to quantitative concepts or prin- 
ciples. 


In filariasis, the statistical distribution of worm 
numbers per host are invariably highly ag- 
gregated, with majority of worms often being har- 
bored by a minority of the host 

population®™. The distribution pattern of parasite 
is similar in both definitive (man) and indefinite 
host (vector) of invertebrate species. The genera- 
tive mechanisms of such patterns are many and 
varied, but include genetic, spatial, and behavioral 
factors. Recently, park (1988) had pointed out 
that the frequency distribution of mf counts in 
samples can be described statistically by a com- 
pound Poisson distribution (the negative binomial 
distribution) which reflects Poisson sampling and 
an underlying (gamma) distribution of mf density 
in the blood of the human population” ’. However, 
it was observed that the distribution of mf in 


human does not fit the negative binomial distribu- 
tion. This discrepancy may be attributed to the 
following reasons: 


i) -Observed negatives are not true negatives as 
proposed by Park et.al due to sampling. 


ii) The observed negatives include, 


a) true negatives who are completely without 
macrofilarial infection and 


b) single sex and prepatent infections. 


Grenfell et al (1989) have refined this analysis, 
proposing the use of a zero truncated negative bi- 
nomial distribution to estimate the considerable 
proportion © of hosts who appear to have no mf in 
their blood. This modification will therefore af- 
fect estimates of the prevalence of infection. Das 
et al have applied this refined analysis and es- 
timated the excess zeros which are due to sam- 
pling bias using the truncated negative binomial 
distribution”. 


2.3.6. Estimation of Excess Zeros in the sample: 


The observed mf frequency distribution for each 
age group is used to fit a zero truncated negative 
binomial distribution by maximum likelihood 
method. The estimated zero probability of this 
distribution, p’(0) can then be used to calculate Z, 
the proportion of the age group who have a true 
zero count (in excess of the proportion due to 
blood sampling), from the equation, 


= [p(0) - p’(0)] / [1 - p’(0)).....(4) 


Where p(0) is the actual proportion of negative 
samples. Finally we can estimate the prevalence 
(R) of infection, corrected for sampling zeros by 
the equation 


R=1-Z 
Example: 
The frequency distribution of mf in 15-19 year age 


class of the people in Pondicherry was fitted to the 
zero truncated negative binomial distribution and 
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wide geographical range’. 


the estimates of parameters » and k were found 
to be 6.26 and 0.27 respectively. The proportion 
of observed zeros (negatives) was 0.895. The es- 
timated proportion of zeros (p’(0)) is the prob- 
ability of observing zero counts using truncated 
negative binomial distribution. These values are 
substituted in equation (4) to get the excess zeros 
‘Z’ due to sampling as 


ie Z = (0.895 -0.426)/(1-0.426) = 0.817 


The corrected prevalence = 1-0.817 = 0.183. 
Similar estimation can be made for all age groups. 


2.4. ESTIMATION OF AGE SPECIFIC 
RATES OF ACQUISITION AND LOSS 
OF INFECTION. 


2.4.1. Introduction: 


In bancroftian filariasis, the age-distribution of in- 
fection is remarkably consistent throughout a 
Typically, both 
prevalence, and intensity rise steadily to a maxi- 
mum in the 15-25 year age class and then decline 
slightly to a stable plateau in adulthood. Such pat- 
terns appear to be ge : egtee by a simple immigra- 
tion-death process” in which the Rate of gain of 
infection (determined by such factors as vector 
biting rate and host resistance) is counteracted by 
the (perhaps determined 
by the host’s ability to terminate infection). Rates 
of acquisition and loss are sensitive parameter not 
only for understanding dynamics of parasite. 
population but also for evaluating control 
programme. Such age prevalence data can be 
analyzed to estimate the rates of acquiring and 
losing the infection by appropriate Caralvtic 
Model developed by Muench * 


2.4.2. Catalytic Models: 


This type of models were originally described by 
Muench” to describe certain chemical reactions 
and later on applied to describe dynamics of in- 
fectious diseases. The basic concepts of Muench’s 
catalytic model is described below. In an endemic 
area for any disease the equilibrium dynamics of 
uninfected (U) and infected (I) human popula- 


a 


tion may be described as follows: 


Mathematically the status of infected population 
(It) over a time period ’t’ can be expressed by the 
following formula 

Ik = Initial infected (Io) + persons 
gained the infection during time 
period’t’ - Persons lost the infection 
during time t 


Changes in the infected population over a time ‘t’ 
can be described as 

I- lo = Persons gained infection - persons 
lost the infection. 


Thus Change in the infected population per unit 
time can be calculated by dividing with total time 
period ’t’. , 


(persons gained infection -persons 
lost the infection)/t 


(It-Io)/t = 


In terms of differential equation it can be ex- 

__ pressed as } : 
dI/ dt = aU -pl (1) 

Similarly the status of uninfected population (Ut) 

after a time period ’t’ will be 


Initial uninfected population (Uo) 
+ persons lost the infection during 
time ’t’ - persons gained the infec- 
tion over a time period ’t’ 


Ur = 


Changes in uninfected population can be 
described as 


persons lost the infection - persons 
gained the infection 


Ut - Uo = 


15 


Thus changes in the uninfected population per 
unit time will be 


(U; - Uo)/t = (persons lost the infection - persons 
gained the infection)/t 


The changes in the uninfected population can be 
expressed in the form of differential equation as 
dU/dt = -BU + ul (2) 
where, 


wis the instantaneous rate of loss of infection from 
an infected individual; 


and 8 is the instantaneous rate of gain of infection 
by an uninfected individual. 


Thus changes in the infection status of a host 
population can be calculated from the rates of loss 
and gain of infection. 


The use of Muench’s models depend upon several 
assumptions. 


a) the most fundamental assumption is that the 
epidemiological situation has remained con- 
stant over the time period t. 


b) the rate of acquiring the infection is same for 
whole population. 


c) the rate of acquisition and loss of infection is 
constant. 


d) the rate of acquisition is the same for those 
persons who have lost the infection as for those 
who never had it or the prior infection does not 
alter the rate at which negative people become 
positive. 


Based on the above assumptions three basic 
models were proposed. 


(a) simple catalytic model,(one way reaction) 


(b) reversible catalytic model,(two way reaction) 


(c) two stage catalytic model (successive reac- 
tion) 


The above three basic models are discussed below 
with examples. 


(a) In the simple catalytic model,the proportion 
of persons showing the sign of infection can 
be expressed by the equation: 


where y is the proportion of persons showing the 
sign of infection after a time period ’t’years in an 
area where the rate of acquisition per person per 
year is B. In this case, it is assumed that the people 
in the endemic area are exposed to the risk of in- 
fection at tie s.me rate, and the method of diag- 
nosis detects all infected persons. 


Eg. Viral infections, where detection of all cases 
is possible. 


This model considers the rate of acquisition and 
not the rate of loss, therefore this model is useful 
in describing such disease where recovery is not at 
all possible like Diabetics and Hypertension. 


(b) Inthe reversible catalytic model, it is assumed 
that the persons in an endemic area are ex- 
posed to a constant risk of infection, conver- 
sion of negatives to positives taking place at a 
constant rate (8), and conversion of positives 
to negatives taking place at a constant rate of 

yw”. Then the proportion of persons (y) show- 
ing infection after a time period of ’t’ is ex- 
pressed by the equation: 


y= {BB+ y)}{1-0P * HI} -....-(4) 


Eg. This model can be applied for diseases like 
malaria and filariasis where conversion of 
negatives to positives and positives to nega- 
tives takes place at a constant rates. 


(c) Inthe two-stage catalytic model, it is assumed 
that persons living in an area endemic for 
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infectious disease are exposed to a constant 
risk of infection, conversion of negatives into 
positives takes place at a constant rate of B, 
positives continue to lose the infection at a 
rate of » and those who lose infections are not 
susceptible to reinfection again. 


Eg. Persons infected with smallpox can become 
uninfected and they will never become in- 
fected again for the same disease. 


In this case, the proportion of persons (y) over a 
time period (t) showing the positive sign of the in- 
fection is expressed by the equation: 


y={B/(B—p)}(eM-e'Pt ) ------(5) 


It is obvious that all of these mathematical models 
are constructed on simple hypotheses neglecting 
the effects of various complicated biological fac- 
tors which differ from disease to disease (like viral 
infection, bacterial, protozoal and helminthic dis- 
eases). Therefore appropriate model need to be 
selected for studying population dynamics of 
pathogen. 


2.43. Application of Catalytic Models in Lym- 
phatic Filariasis: 


Simple catalytic model does not take into account 
the loss of infection, whereas in lymphatic 
filariasis loss and gain of infection is a common 
feature. Secondly, in the case of lymphatic 
filariasis 100% of the infected people can not be 
detected by the present techniques. Hence this. 
model can not be used for describing dynamics of 
lymphatic filariasis. Similarly two stage model as- 
sumes that the reinfection do not occur whereas 
in lymphatic filariasis reinfection is a common fea- 
ture. Therefore this model also can not be used. 
On the other hand reversible catalytic model ap- 
pears to be more appropriate in describing the 
dynamics of W.bancrofti because this model as- 
sumes the chance of reinfection in recovered 
people. Based on the above reasonings the revers- 
ible catalytic model can be used for the estimation 
of the rates of acquisition and loss of W.bancrofti 
infection 


ae ares an ee st 


For age specific rate of acquisition paired sample 


should be collected from an area over a period of 
time. Ideally the area should have reached a 


stable state and acquisition of infection is inter- 
rupted by control measures. The data in the fol- 
lowing example are obtained from a control area 
where transmission is interrupted by efficient vec- 
tor control measures”. 


Example: 
Data source: control area. 


Age Sampled _ Positive No.became 

(yrs) (81&86) in81 negative 

in 5 yrs 
0-5 374 10 6 
6-10 950 55 30 
11-15 796 83 52 
16-20 537 70 44 
21-25 400 56 32 
26-30 426 43 25 
31-40 633 50 36 
41-50 451 35 19 
>=51 310 24 15 
Total 4877 426 259 


___ As the rate of acquisition is nil (i.e. 8 =0) in the 


control area then the equation dI/dt = BU - ul 
will become dI/dt = -yI. By solving the equation 
the rate of loss (uw) can be calculated by the fol- 
lowing formula: | 


yw. = -In(Ir/Io)/t --------- (6) 


where In’ refers to the natural logarithm and It/Io 


refers to the proportion still infected (i.e. from 


1981 to 1986) up to time period ’t’ (in this case t = 
Syears). 


No. of persons + ve both in 1981 & 1986 


No. of + ve persons in 1981 


4 Considering the total cases, out of 426 mf +ve 


persons in 1981 167 remained mf positive in 1986 
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Proportion still infected for 5 years 


= 167/426 
= 0.392 


Instantaneous rate of loss of infection (i.e. x) is 


w = -In(0.392)/5 
= -(-0.9365)/5 
= 0.1873 


2.4.4, Calculation of the life-span of adult worm: 


A more straight forward approach for estimating 
the life-span of adult worm has involved long term 
observation of infected patients no longer ex- 
posed to reinfection either because they have 
moved out of the endemic area or because local 
transmission has been completely stopped. The 
period during which microfilariae were detec- 
table in blood can be calculated by taking the 
reciprocal of (i.e. 1/1) 


The life span of adult worm ( 1/,,) 


= 1/0.1873 
= 5.34 years 


2.4.5. Calculation of gain of infection: 


For estimation of the age specific rate of gain of 
infection paired data can be collected from in- 
dividuals living in an endemic area where no con- 
trol measures are undertaken. Information about 
the status of parasitaemia in each individual has 
to be obtained at two different time points with an 
interval of not less than the prepatent period of 
the parasite under consideration. 


From the following data the gain of infection can 
be estimated as follows: 


The gain of infection = 


No. of persons -ve in 1981 became + ve in 1986 


No. of -ve persons in 1981 


Age(yrs) Sample mf-ve 1981 -ve (’81) to 
+ve (’86) 

0-5 507 494 15 
6-10 1404 1329 42 
11-15 1207 1084 45 
16-20 824 730 38 
21-25 616 535 21 
26-30 671 606 20 
31-40 994 914 24 
41-50 745 685 23 

> 50 557 517 18 
Total 7525 6894 246 


In the Endemic Area, about 0.036 (246/6894) of 
the amicrofilaraemic cohort in 1981 had become 
microfilaria positive by 1986. Now, the change in 
uninfected population over a time ’t’ is given by 
the equation (2) and the solution of this equation 
(ie. proportion of uninfected population over time 
t’) is U/Up, ie. 


Ut/U0 = e-(B + p)t+ {u/(B + )} {1-< — es 
1 ee eee 7 


by subtracting U:/Uo from 1 we can get the 
proportion of persons showing infection after a 
time period ’t’. ie. 

1-Uy/Uo. 


An estimate of the rate of acquisition of infection 
for an equivalent population in which transmis- 
sion continues can then be obtained by substitut- 
ing this estimate of y» in the equation (7). 
However, the real estimate of rate of acquisition 
(B) can be obtained by trial & error method ie. by 
taking the values of B from 0.008 to 2.0 with incre- 
ment 0.0002. For each value of 8, Ut/Uo can be es- 
timated from equation (7). The difference 
between the proportion of observed gain (0.036) 
and estimated gain (1-U;/Uo) can be obtained and 
tabulated in the following manner. 


Fhe value of B at which the sign of difference chan- 
ges first time can be taken as the real rate of ac- 
quisition. 


Estimated B Ut/Uo Difference 
0.008 0.97 0.010017 
0.0082 0.97 0.009386 
0.0112 0.96 -0.00001 
0.0114 0.96 -0.00063 


In practice, however, the apparently uninfected 
and presumed susceptible, population contains a 
proportion (x) of individuals who have a history of 
infection but have become microfilaria negative 
during the progression to disease. These in- 
dividuals are potentially resistant to further infec- 
tion and may introduce bias into the estimation of 
the size of the susceptible populations. This 
potential source of bias may be discounted by as- 
suming that disease rate in amicrofilaraemic in- 
dividuals is constant over 5 years period. 


The proportion of actual susceptible population 
over a time period ’t’ is given by 


St/So = (Ut - xt)/(Uo-xo) 


Where S; is susceptible population at time t- 
So is initial susceptible population 
Ur is uninfected population at time t 
Uo is uninfected population at initial time 
Xt is number of clinical sign cases among 
uninfected population at time t 
xo is number of clinical sign cases among 
uninfected population at initial time. 


The proportion of amicrofilaraemic individuals 
with chronic disease in each age-class should be 
determined by the clinical and blood examination 
of a random, age-stratified sample. 


2.4.6 Calculation of corrected gain: 

Assume that disease pattern is same from 1981 to 
1986. The actual amicrofilaraemic persons (So) in 
1981 is given by 


No. of -ves in 1981 (Uo) - (Uo x proportion of 
disease in -ves (Xo)). 


Data source for corrected rate of acquisition: 


PRPS STS Pe TOF PF eS98eF Fer ee eee csesessecen an2eee 
PPPRI PAS 2 2-22 EEE RENNER 22.5 0c 28 ESeeRenennnedoeeoeesnanatiiinis nnemmee 


—— Cini survey ules 
Age Sample mf-ve -ve (81) Sample prop mf-ve with 
(Yrs ) in ’81 +ve (’86) clinical signs. 
0-5 251 245 9 114 0.00 
6-10 700 658 20 228 0.01 
11-15 587 520 20 216 0.05 
16-20 377 326 22 135 0.11 
21-25 245 208 10 86 0.19 
26-30 225 193 9 56 0.24 
31-50 607 551 21 176 0.33 
> 50 302 277 7 94 0.43 
Total 3294 2978 118 1105 0.14 
considering the total cases mf prevalence is given by 
ie 2978-(2978 x 0.14) = 2561.08 Rate of acquisition x 100 
The actual amicrofilaraemic persons (St) in 5 
years time loss rate of infection 
(i.e. in 1986) is given by This is called point infection prevalence and this 
can be compared with the observed infection 
(no.of mf-ves in ’81- no.with neg. to pos. in 5 prevalence. 


years)-(mf-ves(in ’81) x prop.dis.in mf-ves) 


ie. (2978-118)-(2978 x 0.14) = 2443.08 2.5. ESTIMATION OF RISK OF DEVELOP- 
ING DISEASE IN RELATION TO 
So, the corrected gain is given by MICROFILARAEMIA STATUS, AND, 
TRUE INFECTION PREVALENCE IN 

1 - (actual amicrofilaraemic in ’86(St) A POPULATION. 
In lymphatic filariasis it was observed that most of 
actual amicrofilaraemic in ’81(So)) the people with chronic manifestations are 
microfilaria (mF) negative on blood examination. 
ie 1-(2443.08/2561.08) = 0.04607 Since individuals who had ever been infected are 


the potential source of developing the disease, es- 
For males the loss rate of infection (m) is found to tablishing a relationship if any, between the infec- 


be 0.17746. Having the corrected gain and loss tion status and the consequent progression of 
__ rate of infection for males the corrected rate of ac- disease in an infected individual, may provide an 
quisition can be obtained by iteration using equa- estimate of the risk of an infected individual 
tion (7) as done for the estimation of B developing the clinical manifestations . The main 
(uncorrected for disease). The corrected rate of aim of this analysis is to estimate the risk of acquir- 


acquisition is found to be 0.0144. The estimated ing disease in persons who had been 
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2.5.1. Estimation of the proportion ever infected 
using catalytic model: 


In the reversible catalytic model, 
y=[B(B+p)J[1-e(— BrMM |. (2) 


assuming that there is no loss of infection i.e., » = 
0, the proportion of population that would have 
remained mf positive for ever without losing the 
infection, (ct) over a period of time ‘t’ is given by 


where ‘’ is the instantaneous rate of acquisition 
of infection. Since the gain rate is age dependent 
(see section - 2.4), the cumulative infection 
prevalence has to be estimated for each age class 
independently. This can be obtained by correct- 
ing the proportion who had infected in the pre- 
vious years of their life. Thus, if utn is the 
proportion uninfected at time tn, then the propor- 
tion ever infected at time tn + 1 is equal to, 


Ctn+1= 1-urne Pin + 1" + 1'y a ) 
Where, utn = 1-ctn. 
Ctn = proportion ever been mf positive at time ‘tn’ 


B tn = rate of acquisition of infection between the 
time periods ‘tn-1’ and ‘tn’ 


Here the two time periods are the mid points of 
the age classes. 


2.5.2. Estimation of the risk of developing clini- 
cal symptoms: 


Using the age specific rates of loss and acquisition 
of infection (‘yr’and‘f;’), the infection prevalence 
(corrected for disease) for age ‘t’ is calculated as 
follows: 


Pt = Bt/pt «ee. (4) 


The proportion of the people who are at the risk 
of developing the clinical symptoms of filariasis 
(Rt) can be obtained by subtracting the corrected 
prevalence from ct. 


1.€. Rt = Ct ~ Pt cvveeceeeeee (5) 


This proportion ‘R;’ of people who had been 
microfilaria positive earlier and become ap- 
parently negative now, (who appear to have 
recovered from infection) is compared with the 
chronic disease prevalence. It was seen that in 
Pondicherry particularly in the male population, 
the recovery (risk) rate resembled closely the 
chronic disease prevalence suggesting the pos- 
sibility that the recovery process triggers off in- 
flammatory reaction and disease. 


This corroborates with the immunological theory 
of lymphatic filariasis which says the response of 
the immune system to the infection results in the 
clearance of mF and death of adult worms there- 
by creating blockades in the lymphatic system 
which in turn leads to clinical manifestations”. 
The inference made from this analysis is that the 
people who have become mF negative from being 
positive are more at risk of developing disease 
than the ones who are microfilaraemic. 


Example : 


The reinfection data was used to estimate the age 
specific rates of acquisition (8) and loss (1) of in- 
fection. The methods of estimation of these two 
parameters are explained in the previous chapter. 


The age specific rates of acquisition and loss of in- 


fection for males in the first few age classes are 
given below: 


Age class Rate of Rate of 
(in years) Loss Acquisition 
0-5 0.0575 0.0087 
6-10 0.1685 0.0091 
11-15 0.1788 0.0120 


The proportion who have ever been infected in the 
first age class 0-5, can obtained using equation (3) 
as follows: 


i pt i nt i i 


———————— i 


here, to = 0, tl = 2.5, Bt: = 0.0087 and Uto=1 
025 = 1~¢ (-0.0087x25) 
= 0.0215 
For 6-10 age class, 


here, tl =2.5, t2 = 8.0, Bi2=0.0091 
and uti = 1-0.0215 


cg = 1-[0.9785 x ¢ 0091 x 8:25) 
= 1-0.9307 
0.0693 


Proceeding this way, we can calculate the propor- 
tion ever positive for all the age classes. 


From the proportion ever positive for mF, we next 
estimate the proportion of the people who are at 
the risk of developing the clinical symptoms of 
filariasis. Using the age specific rates of loss and 
gain of infection (‘Br and ‘yt’), the prevalence of 
infection for the age class 6-10 can be obtained 
using equation (4). 


ie. ps = 0.0091/0.1685 

= 0.0540 
The proportion of ever infected people who are at 
the risk of developing the clinical symptoms of 
filariasis (Rt) can be calculated using equation (5). 


i.e. Rg = 0.0693 - 0.0540 
= 0.0153 


Similarly the estimate of risk of developing the 
clinical symptoms of filariasis can be obtained for 
all age classes. 


2.5.3. Estimation of true infection prevalence 
using the probability of acquiring a male or 
female worm: 


The transmission dynamics of filariasis involves 
three factors, namely, (i) the human host, (11) the 
parasite and (iii) the vector. The human popula- 
tion in an endemic area can be either infected or 
uninfected. The infected population are further 
subdivided into the following: 
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a) microfilaraemic, 

b) individuals with very low worm burden, 
c) single sex (worm) infections 

d) individuals at prepatent period and 

e) diseased 


For any disease, one of the measures of transmis- 
sion is its infection level in the population. Infec- 
tion level is defined as the proportion of persons 
harboring the parasite causing the disease in a 
population. Conventional method of estimating 
the disease\infection prevalence identifies only 
the microfilaraemic population and ignores all 
other groups (from (c) to (e)) due to following 
technical difficulties. 


In Lymphatic filariasis, a person is said to be in- 
fected the moment the infective larvae enters his 
body. The conventional blood sampling, can not 
detect the infected persons during the prepatent 
period. It is also less sensitive to detect carriers 
with ultra low level mf in their blood. Although 
membrane filter concentration technique 
facilitated the detection of low density carriers of 
microfilariae, its application in large scale survey 
is impracticable due to various factors*. Since the 
parasite is of dioecious in nature, these techniques 
can not detect persons with single sex (worm) in- 
fection. Therefore, it is practically not possible to 
estimate all the ‘ infected’ persons in a population 
for the above mentioned reasons. However, an es- 
timation of infection prevalence which takes ac- 
count of few of these factors will be useful to 
define the public health importance of filariasis in 
an area. 


In section 2.3, it has been shown that the fitting of 
truncated negative binomial model to the fre- 
quency distribution of various mf counts in the 
population also accommodates the single sex and 
single worm infections while correcting the ob- 
served prevalence of infection”. In this analysis 
the true infection prevalence (including single sex 
(worm) infections) in the population is estimated 
using the probabilities of occurrence of male or 
female worms based on some assumptions. The 


method uses the age specific rates of loss and gain 
of infection calculated using the reversible 
catalytic model (described in detail in section- 
2.4). Consequently this analysis gives an account 
of the changes in the mF prevalence due to single 
sex (worm) and its significance in disease develop- 
meiit (as shown by Hairston and Jachowski in the 
case of aperiodic bancroftian filariasis”). The 
basic assumptions for calculation of probability of 
acquiring a male or female worm are: 


(i) | Theevent of acquisition of a male worm by 
an individual is independent of him\her 
acquiring a female worm. 


(ii) | Any infected individual will harbor equal 
number of male and female worms i.e., 
he/she will harbor equal ratio of male and 
female worms. 


Violation of these assumptions might lead to er- 
roneous results and there by making the inter- 
pretations difficult. 


Let ‘p’ be the proportion of persons positive for 
infection as obtained from the microfilaria survey. 
This proportion of people harbor both male and 
female worms which is evident from the presence 
of mf. Based on the assumption (i),the proportion 
positive for infection should be equal to the 
product of the proportion harboring male worm 
(pm) and proportion harboring female worm 
(pr).(Since pm and pr are independent). i.c., 


P = PmXx Pf 
By assumption (ii), 
Pm = Pf=Pc 


Therefore, , 
ae 
where pc is the constant proportion containing 


either male worms or female worms. Therefore it 
follows that 


Pm = pf = sqrt(p) = pe 


Based on the same arguments, substitute the 


(ct) in the place of 


proportion ever positive for mf 
proportion positive (p). We then calculate the 


proportion who were ever infected with either 
male worm (ym) or female worms (yf) by taking 
the square root of the proportion who had ever 
been positive for mf over a period of time ‘ Le., 


Ym = Yf = SQTU(Ct).--eeereeeeeee(6) 


Then the proportion of persons never infected 
with female worms is given by 


(1-yr) 


Since there is an equal sex ratio among the worms, 
the proportion which escaped the infections due 
to male worms also is given by 


where yp is the proportion which escaped all the 
infections (due to both male and female worms) 
i.e., the one which is uninfected. The proportion 
ever infected with worms of either sex, ‘y;’ can be 
estimated by the following expression 


Now this ‘true’ infection prevalence can be calcu- 
lated for different age classes and is compared 
with the disease prevalence. It was found that in 
Pondicherry there was no resemblance between 
the two suggesting the negligible role of single sex 
(worm) infections in disease development in the 
case of nocturnally periodic bancroftian filariasis. 
However interpretations based on the estimates 
of true infection should be done with enough care 
since all these calculations are based on the fact 
that the assumptions (i) and (ii) are true under all 
situations. 


It should be also noted that even this estimate of 


true infection prevalence does not include the 


persons in prepatent period. Therefore to include | 


those individuals also, it is necessary to measure 
the number of L3 larvae that successfully enter the 
human system thorough immunological methods. 


Example : The proportion ever infected (cr) is 
0.0215 for the age class 0-5. The true infection 


prevalence can be estimated as follows: 


Using equation (6), the proportion ever infected 
with female worms, mated or not is given by 


yf = sqrt(0.0215) 
= 0.1467 


Then the proportion of persons never infected 
with female worms is given by, 


= 1-0.1467 
= 0.8533 


By equation (7), the proportion which escaped 
the infections of both male and female worms is 
obtained as below 


= (0.8533) 
= 0.7281 


Then the true infection prevalence (persons ever 
infected with worms of either sex, mated or not) 
is estimated using equation (8), 


= 1-0.7281 
= 0.2719 


Proceeding in a similar manner, the true infection 
prevalence can be estimated for all age classes. 


2.6. SPATIAL DISTRIBUTION OF 
PARASITE CARRIERS 


In section-2.3 the negative binomial probability 
distribution was used to describe the distrivution 
of worm aggregation in the human host. Since 
clustering of infected individuals in households is 
potentially an important consideration, it is neces- 
sary to study the spatial distribution of mf carricrs 
in relation to size of households. Several 
published reports are available in this aspect 437" 

? Ottesen et al.. .,(1981) have also reported that 
the familial clustering is also due to genetic in- 
fluences on the susceptibility to bancroftian 
filariasis’. Kimura et al ., demonstrated the 
clustering of mf carriers by considering only those 
who have been examined (irrespective of family 
size) for microfilaraemia. Since in their study, all 


____ the members were not examined, the proportion 


infected in each household is biased to the num- 
ber examined’*. Vanamail et a/ . has extended 
this analysis for each household size, thus any bias 
resulting from differing proportions of infected 
individuals in household of different sizes was 
avoided*’. In this section the method of collect- 
ing and analysing data on the spatial distribution 
of mf carriers are discussed, 


METHODS:- 


In an endemic area the house wise enumeration 
(age & sex) of persons should be carried out prior 
to any parasitological survey. During enumeration 
time it is also important to identify the various fac- 
tors that are supposed to influence the transmis- 
sion dynamics of the disease with respect to each 
house. 


Enumeration Form for Microfilaria 
Survey 


Format: 


House Age/ Breeding History of Other 
No. Sex Source cases factors 


After completing the enumeration with all 
relevant informations, the parasitological survey 
should be carried out. During this study as much 
as possible all persons in the target house should 
be covered. Equal weightage should be given for 
all size of households. 


Tabulation: From the enumeration and 
parasitological survey data the fol- 
lowing Table should be completed. 


Family No. of houses 

Size Enumerated Full y Partly 
covered covered 

1 

2 (* Example) 

3 

4 

i 250 117 70 


2.6.1. Analysis: 


For clustering analysis it is more reliable to con- 
sider the data on full participated houses in the 
study. The results of parasitological survey from 
the full participated houses can be obtained for 
each family size as follows. 


Family Houses with +ves 
Size 0 1 2 3 4 
1 

2 

3 

4 

_ a2 30 4 1 0 
6 

7 

8 

9 

* - example 


For each family size the frequency distribution of 
parasite carriers can be obtained. 

Example: Frequency distribution of parasite 
carriers for family size 5. 


a a ee ee ee 


In order to see the distribution Pattern of infected 
persons in relation to family size, statistical dis- 
tributions can be fitted to the observed data””. If 
the distribution of infected individuals is random 
then Poisson probability model will describe the 
observed trend in a significant manner. The form 
of Poisson distribution is 


Pa) =e¢ px! x = 0,1,2....... 


Where, ’w’ is the mean of Poisson distribution 


If the Poisson distribution do not describe the 
trend adequately then the distribution of infected 
persons deviate from randomness and it implies 
that the distribution of infected individuals is ag- 
gregated or clumped. Since negative binomial dis- 
tribution is one of the clumping distributions, 
aggregation phenomenon can be tested by fitting 
negative binomial probability statistics. The 
method of estimation of parameters of this dis- 
tribution are described in section 2.3. The testing 
of randomness or clustering can be carried out for 
each family size. 


Fitting of Poisson distribution to 
the above frequency data for fami- 
ly size 5. 


Example: 


0 82 82.4 0.00207 
1 30 28.9 0.04344 
2 4 5,3 0.08596 
3 1 0.6 

Total 117 117.0 0.13147 


Here calculated Chi-square value is 0.13147. 
While comparing this with the tabulated value for 
1 degree of freedom (i.e number of observations _ 
in the Chi-square column-1-number of 
parameters to be estimated) the level of goodness 
is more than 0.05 (ie P >0.05). The goodness of 
fit indicates that the distribution of infected in- 
dividuals is random. In certain cases the family 
sizes can be added into 3 or more groups due to 
less number of observations in each family size. 
Adopting the above two statistical approaches 
one can get two types of households i.e. the 
households under the family size at which ran- 
domness appear and the households under the 
family size at which clustering phenomenon oc- 
curs. The significant causes for difference in dis- 
tribution pattern can be analyzed further by 
considering age and gender distribution of the 
population in the two types of households and also 
other possible factors that influence the transmis- 
sion dynamics of the disease. 


3. ENTOMOLOGICAL DATA 


3.1. Introduction: 


Transmission dynamics of any disease involves the 
study of interaction between the parasite with man 
and the vector. In the previous chapters we have 
described the methods of sampling the human 
population for studying the distribution of 
parasite in the host and their practical implica- 
tions in estimating the prevalence and the estima- 
tion of fecundic life span and age specific 
incidence of the parasite in the human host. Fur- 
ther the relationship of microfilaria infection with 
disease prevalence are explored. In this chapter 
the methods of sampling the larval and adult stage 
of the vector population (the intermediate host) 
and the parasite distribution in the vector, and its 
relationship with parasite distribution in human 
population are described. 


The various methods of sampling and estimation 
of the larval and adult population are well 
described in many text books’ *"’’”. A recent 
WHO publication has well documented the 
various methods of sampling and estimating the 
adult population of C.quinquefasciatus 2 We 
begin this chapter with the methods of sampling 
the intermediate host of the parasite. The life 
cycle of the intermediate host ( C. quinquefas- 
ciatus ) involves two stages, namely, the immature 
stage and adult stage. Since the immature and 
adult stages are living in entirely different habitats, 
the sampling methods also varies accordingly. 


3.2. SAMPLING THE LARVAL POPULA- 
TION: 


The purpose of larval sampling is of many fold. 
The following are the few among many objectives: 


1. To monitor and give feed back to control 
operations. 


2. Tostudy the dynamics of the larval population. 


3. To study the relative importance of species 


25 


— in the transmission of the dis- 
ease 


4. To study the relative importance of different 
types of habitats for breeding of a particular 
species ~. 


The choice between one time survey or con- 
tinuous surveillance of the immature population 
will depend on the objectives”. If the objective is 
to give feed back to larval control operations then 
a continuous surveillance is necessary. The 
second and third objectives may be fulfilled by one 
time or continuous surveillance of the larval and 
adult population. An account of the need for a 
sampling design and the advantage of using 
stratified random sampling are described with 
suitable examples in " A note on Sampling 
methods " by S.J. Poti“ . Though the methods we 
describe here are specie to the larval sampling of 
C. quinquefasciatus breeding habitats in an urban 
situation, it may be applicable to other species of 
mosquitoes whose breeding habitats are similar in 
nature. 


3.2.1. The sampling design: 


The major breeding sources of C. quinquefas- 
ciatus are cess pits, cess pools, unused. wells and 
drains in an urban situation’®*?, The study area 
can be broadly classified into low, moderate and 
heavy breeding potential areas. Each of the above 
strata is comprised of substrata like villages. All 
the villages can be serially numbered and the 
habitats identified by house numbers or by giving 
some identification numbers to each type of 
habitat. For this purpose an initial survey may be 
carried out to enumerate the cess pits, cess pools, 
wells and drains. The form-1 may be completed 
for each habitat in each village. Since the breed- 
ing habitats may be widely scattered it is impos- 
sible to cover the entire area in each of these 
stratum. Alternatively a sample of habitats may 
be randomly selected in each stratum and dipper 


or bucket samples are taken only from these 
habitats. For each type of habitat the form - 2 may 
be completed. An estimate of the mean of the 
population per fixed number of dips for each type 
of habitat can be calculated by the following for- 
mula: 


m=XNjyui/2Ni Where, 
mj = sample mean of ith stratum = 
> xi/d Di 
Xij = Total number of larvae per fixed 


number of dips from jth unit (cess 
pits or wells or cess pools) of ith 


stratum. 

ni = Total number of units sampled in 
the ith stratum. 

ha Total number of units in the ith 
stratum. 


The variance of the mean ‘m’ is given by Cochran 
as, 


variance of (m) = £(Wi Si/nj)’ 
where, Wj = Ni Ni 

Si? = [1/(ni-1)][2 (xj - mi)”] 

= Sample variance of ith stratum. 
A comparison of different habitats can be made 
by taking the weighted average of estimated mean 
population in each habitat with weight being equal 
to the proportion of area actually breeding to total 


area of that habitat in each stratum”. 


FORM - 1.: Breeding habitats enumeration form. 


Type of Habitat Date: 
S.No. Village Total No..of Surface 
code habitats area (Sq.m) 
5 10 120 500 


FORM -2.: Proforma for Larval Sampling 


Village Code: Habitat: Date: 


0 OO 66 OOO SS SSS 66S OS 8S 888888 BEE OOS 8 Ee 


Habitat No.of. No. of larvae in stage Total 


LD.No. dips. I IM Il IV P Larvae 
taken 
15 6 48. 2 22: 2 eee 


3.3. SAMPLING THE ADULT POPULA- 
TION: 


As it is not possible to examine the entire 
mosquito population in a large area, a repre- 
sentative sample has to be collected over a fixed 
period of time. The purpose of sampling the adult 
mosquitoes can be broadly classified as follows: 


1. To study the population dynamics of the vec- 
tor. 7 


2. To study the transmission dynamics of the 
parasite from vector to man. 


3. To give feed back and evaluate the vector-con- 
trol programmes. 


Adult mosquitoes can be sampled from the rest- 
ing or biting population. Collection of resting or 
biting mosquitoes by hand catch provide more — 
representative samples of the population. Rest-. 
ing or biting mosquitoes are collected by hand 
using either aspirators or a test tube for a fixed 
time period from the same houses or baits. The 
data obtained can be expressed as the average 
number of mosquitoes collected per man per hour 
as an index for comparing the seasonal fluctuation 
in the population. Though sampling the biting 
population would be ideal from epidemiological 
point of view, it is less practicable due to the fol- 
lowing constraints: 


1. Getting baits for mosquito collection. 


oy 


2. Ifthe population is of periodic in nature, then 
a biting sample should coincide with the peak 
biting hours. Since C. quinquefasciatus are 
known to be nocturnally periodic, a sample ob- 
tained from early or late hours may not provide 
a representative sample. Hence a whole night 
biting collection has to be made which is prac- 
tically difficult. 


Since the resting mosquitoes include a random 
mixture of fed as well as unfed mosquitoes and it 
is practically possible to collect mosquitoes during 
day time, it would not be biased if estimates are 
obtained from a sample of resting population. 
This could be true at least, if not all, for mosquito 
species with nocturnally biting rhythm. Further a 
whole night biting collection is useful only for 
determining or identifying the local vector species 
and their periodicity in nature. 


3.3.1. Estimation of entomological Parameters: 


The entomological parameters commonly used in 
the evaluation of control programmes are, the 
relative biting or resting densities, survival of the 
vector population, infection and infectivity rates, 
the annual infective biting rate (AIBR) and the 
annual or monthly transmission potential (ATP or 
MTP) of the vector population. WHO (1984) in 
its fourth report on lymphatic filariasis stated that 
"a comparison of these parameters, especially the 
AIBR and ATP, before and after control opera- 
tion is valid, provided the sampling and dissection 
methods are standardized,". Some of these impor- 
tant parameters are described» below. 


Annual Biting Rate (ABR): 


This is the number of bites received by one man in 
one year. This parameter is estimated from 
direct 12 hours human bait landing catches car- 
ried out once in a fortnight or a month for a year. 


Annual Infective Biting Rate (AIBR): 


This is expressed as the number of infective 
mosquitoes biting a man in one year. The product 
of the infective rate in each catch and the es- 
timated fortnightly or monthly biting rate gives the 
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infective biting rate for that period. The AIBR is 
obtained by the summation of this parameter for 
a year (24 fortnights or 12 months accordingly). 


Annual Transmission Potential (ATP): 


It is the number of infective larvae that could have 
been received by one person in a year. This 
parameter is calculated using the following for- 
mula. 


ATP = aXbXc where, 
a= Estimated number of mosquitoes biting a 


man in one year (This is calculated as per 
man hour biting density X 12 hours X 365 


days). 

b = _ Proportion of mosquitoes infective (from 
biting collection only). 

c= Number of infective larvae per infective 


mosquito. 


Further, we describe some new methods for es- 
timating the survival rate of adult mosquitoes and 
the relevance of currently available parameters 
for measuring the level of transmission. 


3.3.1.1. The survival rate of vector: 


The survival rate of vectors is an important com- 
ponent in determining the efficiency of transmis- 
sion. This is the basic for determining the 
expected infective life span of the vector. One of 
the commonly used methods for estimating sur- 
vival rate is based on the mean purous rate. Field 
collected mosquitoes are dissected and graded 
according to their parity status as parous or nul- 
liparous etc. Then the survival rate (S) is calcu- 
lated using the following formula: 


S = nth root of proportion parous (P) 
where 
Total parous 


P = Parous rate = --------------------- 
No. dissected 


n =the interval between two blood meals. 


The standard error may be estimated by the fol- 
lowing formula: 


Standard error (S) = sqrt [S (1-S)/no.dissected] 


The above formula makes the following assump- 
tions: 


1, The recruitment rate of nulliparous 
mosquitoes is constant between two time 
points of survey. 


2. The population is stable., and 


3. The duration between blood meals is constant 
for all stages of the mosquito. 


However the recruitment rates and interval be- 
tween blood meals are seldom constant in nature 
and consequently the changes in the proportion 
parous may bias the estimate of the survival rate. 
Based on the above reasonings new technique has 
been developed for populations which are rapid- 
ly changing over time*°. These authors have ex- 
pressed the dynamics of the biting mosquito 
population by two equations as below: 


Mt = PuTt-u 
Te = Net+Mt 


where T; is the total biting population on day t, 
composed of N; nulliparous and Mt parous in- 
dividuals. The proportion Py is the survival rate 
per gonotrophic cycle, which has a mean duration 
of u days. Then an estimate of the Py is obtained 
by the following formula, 


Py = 2 Mri/2Tii-u 


Where Miis the parous sample in the total sample 
Fti collected on day tj, and there are m pairs of 
sampling occasions separated by u days. This 
method assumes that the time series M(t) and T(t- 
j), lagged by j days are correlated. Cross correla- 
tion coefficients are estimated with varying time 
lags and the lag at which the peak correlation oc- 
curs is the best estimate of the mean length of the 


cycle. The introduction of the time lag u is as- 
sociated with the gonotrophic cycle of the biting 
population. This time lag solves the problem of 
varying recruitment rate, which biases the mean 
parous rate formula. However the integer u, the 
mean length of the cycle, is unknown and it, too 
must be estimated from the data. 


3.3.1.2. Distribution of parasite in the vector 
mosquitoes.: 


The distribution of microfilaria in the human 
population has been adequately described either 
by the negative binomial or lognormal probability 
functions’”’. However, the fitting of lognor- 
mal probability distribution to the observed 
microfilaria distribution did not have any biologi- 
cal explanation. Further modelling a discrete 
data (microfilaria counts) by a continuous prob- 
ability distribution is not reasonable. Grenfell et 
al. (1990) and Das et al.”” have explained the 
biological mechanisms that generate the over dis- 
persion of parasite in the human population and 
the consequence of describing this by a negative 
binomial distribution. : 


Das etal. (1991 yet have recently observed that the 
distribution of parasite infection by vectors of lym- 
phatic filariasis (as measured by mf counts) can be 
described by the same statistical model as that 
proposed for blood sampling mf by Grenfell et 
al. ( 1990). In particular, the mf count distribu- 
tion in nulliparous mosquitoes is remarkably well 
described by the (zero-truncated) mixed Poisson 
distribution, which was used to fit the equivalent 
pattern of mf counts in humans. 


The similarity of the observed vector and human 
distributions, therefore indicates either that 
counts from mosquito dissections reflect a higher 
sampling efficiency, or that there may be some 
concentration of parasites by the vector. The 
overdispersion analysis also confirms that parasite 
distributions become less overdispersed with 
parasite stage, reflecting a lower proportion of 
vectors with high parasite loads. This provides 
evidence of density dependence in the host- 
parasite relationship, either in parasite survival or 
parasite-induced host mortality. 


4. CONCLUDING REMARKS 


A series of important papers were published on 
reviewing quantitative techniques that are often 
used for analyzing data from filariasis control 
programmes” 0” The authors discussed and 
criticized the methods used for estimating the 
microfilaria density and periodicity until the end 
of 1960 and provided some other methods. 
However, none of the workers seem to have 
pointed out the bias while estimating parameters, 
which is the function of sampling techniques 
among others. Aikat ef g/ emphasized the need 
for probability sampling and provided a stratified 
two-stage probability sampling design” 1 But 
their design did not take into account the factors 
that influence the estimation of parameters. 
However, from their study valid inferences were 
drawn to make use of it in the subsequent surveys. 
The sampling design proposed here takes care of 
practical constraints the authors have faced while 
evaluating a control programme in Pondicherry. 
The design was prepared by giving due considera- 
tion to the factors that would affect the estimation 
of parameters while evaluating any control 
programme. 


The assessment of filariasis control programmes, 
either by chemotherapy or vector control, is 


29 


based on clinical, parasitological, and 
entomological methods!. The parasitological as- 
sessment is mainly based on the microfilaria 
prevalence and intensity ina community. The de- 
gree of reduction is an indication of effectiveness 
of the control measures. This is particularly so in 
drug control programmes, which directly affect 
the fecundity of adult worms. However, the assess- 
ment based on this parameter in areas where only 
vector control measures were carried out, leads to 
paradoxical situations. This could be avoided by 
applying more sophisticated statistical techni- 
ques. 


Though the understanding of parasite distribution 
in human host and its relationship with disease are 
well established, the knowledge on parasite dis- 
tribution in the vector and its transmission 
dynamics to man and vice versa is at infant stage. 
The quantitative relationship between infection 
rates and intensities in the human population and 
the amount of transmission of the parasite to the 
development of microfilaria and of disease are es- 
sential for planning and monitoring of vector con- 
trol programmes”. Undoubtedly these are the 
unresolved puzzles which require further im- 
munological and statistical/mathematical techni- 
ques. 
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